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ABSTRACT 

The genus of the iso-density contours is a robust measure of the topology of large scale structure, and 
it is relatively insensitive to nonlinear gravitational evolution, galaxy bias and redshift-space distortion. 
We show that the growth of density fluctuations is scale-dependent even in the linear regime in some 
modified gravity theories, which opens a new possibility of testing the theories observationally. We 
propose to use the genus of the iso-density contours, an intrinsic measure of the topology of large-scale 
structure, as a statistic to be used in such tests. In Einstein's general theory of relativity, density 
fluctuations are growing at the same rate on all scales in the linear regime, and the genus per comoving 
volume is almost conserved as structures are growing homologously, so we expect that the genus- 
smoothing scale relation is basically time-independent. However, in some modified gravity models 
where structures grow with different rates on different scales, the genus-smoothing scale relation should 
change over time. This can be used to test the gravity models with large scale structure observations. 
We studied the case of the f{TZ) theory, DGP braneworld theory as well as the parameterized post- 
Friedmann (PPF) models. We also forecast how the modified gravity models can be constrained with 
optical/IR or redshifted 21cm radio surveys in the near future. 



1. INTRODUCTION 

The large scale structure has long been a major 
source of information on cosmic evolution. The most 
frequently used statistics of the large scale structure, 
such as the two-point correlation function and the den- 
sity power spectrum, and especially the baryon acous- 
tic oscillation (BAO) signatures in the power spec- 
trum have been widely used in tests of cosmological 
models and precisional measurements of cosmo l ogical 
parameters (iTegmark et all 120061 : iPercival et~all 120071 : 
iGaztanaga et al.ll2009D . However, these two-point statis- 
tics do not exhaust all of the information content about 
the large scale structure-at least not in the non-gaussian 
case. There are also some limitations on the application 
of these statistics. As the perturbations grow larger, 
their evolution becomes nonlinear, resulting in a dis- 
torted power spectrum and the emergence of gravity- 
induced non-gaussianity. This has a crucial influence 
in the extraction of the absolute BAO scale: at fc = 
0.1 ~ 0.15Mpc/ft. where the flrst BAO peak is located, 
the non-linear evolution of the power spectrum amounts 
to 2% to 5% diff erences between z = and z = 0.5 
(|Kim et. al.l 120091) . Accurately modeling the non- linear 
evolution is highly non-trivial, and practical application 
in actual model test/parameter determination also re- 
quires fast and easy-to-implement algorithm. Although 
progresses have be en made in more sophisticated pertur- 
b ative expansion (iCrocce fc Scoccimarro 2006a, b, [2003; 

fatsubara"2008a"b'; 'Smith et al . ,2007i) or re construction 
admanabhan et al. 2009; Noh et al]|2009[ ) procedures 
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to correct for the effect of non-linear evolution, a precise 
yet practical method to take into account of these effects 
has yet to be developed. Other systematic effects, such 
as the clustering bias and the redshift distortion could 
also affect the outcome of the measurement. Such effects 
may significantly limit the accuracy of cosmological tests 
for future redshift surveys, where the statistical error on 
the BAO scale would be at the one percent level. 

An alternative approach to the correlation function 
and power spectrum is the topology of the iso-density 
contours, in particular the genus of these contours. This 
offers another way to characterize the statistical prop- 
erty of the large-scale structure, and is insensitive to 
the systematic effects discussed above, since the intrin- 
sic topology does not change as the structures grow, at 
least not until th e iso-density c ontour s eventually break 
at shell crossing ljPark fc Kim I [20T0I : iPark et all [2005). 
This provides a robust statistic for cosmology. Com- 
pared with the power spectrum, the genus chang es only 
by 0.5% between z = 0.5 and z = (Park1 l2011[ ). This 
is particularly true when using the volume fraction for 
the threshold levels to identify the contou rs. Accord- 
ing to the second order perturbation theory (jMatsubara I 
Il994f ). there is no change in the genus at median density 
threshold due to the weakly nonlinear evolution, because 
the iso-density contours enclosing a given fraction of vol- 
ume do not change as long as the gravitational evolution 
conserves the rank ordering of the density. Similarly, a 
monotonic clustering bias does not result in any differ- 
ence in the iso-density contours, and a continuous co- 
ordinate mapping into redshift space does not affect the 
genus statistics either. In ACDM model with general rel- 
ativity, the genus is conserved during the different epochs 
in the linear regime of evolution. Therefore, by observing 
the genus curve at the different redshifts and smoothing 
scales, one can use it as a robust standard ru ler for cos- 
mological measurements ([Park fc Kim |[2010l ). 

In the present paper, we consider the phenomenologi- 
cal consequences of alternative gravity theories. A num- 
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ber of such models, e.g. the /(TZ) theory (iCarroll et al.l 
200l ICapozziello et al.l [200a INoiiri fc Odintsov 1 120031: 
Song et aLr2007aD and the DGP brane- world model 



(jDvali et al. I2000t ). have been proposed to explain the ac- 
celerated expansion of the Universe. Unlike the general 
relativity, where the growth of the density fluctuation is 
at the same rate on all scales during the linear evolution, 
the modified large-scale gravitational forces induced by 
extra scalar field can in general introduce new scale de- 
pendence in the growth of structure, and therefore distort 
the genus curve in a time-dependent way. This provides 
us a new tool to distinguish GR from its various alterna- 
tives. It probes the combined effects of background ex- 
pansion and the growth of the structures, and therefore 
is able to break the degeneracy between the dark energy 
and modified gravity models from expansion data alone 
(jLinder 1 120051 : 'Zhang et al. 20071. Furthermore, such a 
measurement also has the advantage that it is insensitive 
to nonlinear systematics. This is particularly useful be- 
cause for some modified g ravity models, such as the on es 
invoking the chameleon (jKhourv fc Weltman I l2004allHl 
mechanis m or the Vainshte in mechanisms (jVainshtein I 
[1972; Dcf favet et al. I l20?)l . it is non-trivial to treat 
the effect of nonlinear evolution on the power spectrum, 
and a conservative treatment on nonlinearity may only 
yield fairly weak c onstraints on gravi ty models with even 
strong signatures (ISong et al.ll2"007bh . The genus, on the 
other hand, is almost conserved during much of the evo- 
lution, thus allowing us to avoid dealing with such prob- 
lems in the nonlinear regime. Instead, we may make the 
calculation in the linear regime, and expect it to be al- 
most unchanged until very late into the nonlinear regime. 

Below in Sec. II we briefly review the calculation of 
the genus density and the cosmological perturbation the- 
ory in the parametrised post-Friedmann framework. In 
Sec. Ill we discuss how the genus density changes in 
the various modified gravity models. We forecast the 
prospects of observation with the future redshift surveys 
in Sec. IV and conclude in Sec. V. An estimate of sta- 
tistical uncertainty of the genus measurement is given in 
the Appendix. 

2. THEORY 

The genus density G is given by the genus g, i.e. (the 
number of holes - number of isolated regions) divided 
by the volume V. Smoothing the tracer distribution 
by a Gaussian filter with scale Rg, then the iso-density 
contours of smoothed field can be identified for a given 
variance- normalized density threshold v = S/a{Rg) or 
volume fraction. The genus curve, i.e. the genus per 
unit volume of the iso-density contours as a function of 
such threshold, could be measured. 

In the particular case of Gaussian fluctuations, the 
genus curve ca n act ually be calcu lated analytically 
(jHamilton et al.lll986l: fDoroshkevich II197O0 : 



(1) 



and the amplitude A is related to the power spectrum 

by 

1 / Jd'kk'Pik)WikRg) Y/' 

in^ 33/^\ Jd^k P{k)W{kRg) ) ' ^' 



where W is the window function for smoothing. The 
threshold level v is related to the volume fraction / on 
the high-density side of the density contour surface via 



/ = 



(3) 



Although the genus is related to the power spectrum 
in the Gaussian random field case, its measurement is 
completely independent of the power spectrum mea- 
suremen t, but rather by me ans of the Gauss-Bonnet 
theorem (|Hamilton et al.lll986D 
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where g is related to the integration of the Gauss cur- 
vature K = l/(rir2) over the surface, ri and r2 are the 
principal radii of curvature at the integration point. In 
practice, after smoothing the galaxies samples or simu- 
lation data into continuous density field, one first makes 
a fine tessellation of the whole volume of the space with 
small polyhedrons, e.g. cubes or truncated octahedrons, 
the iso-density surfaces can then be approximated by the 
polyhedral surfaces. Since the only nonzero contributions 
to the integrated curvature of the polyhedral surfaces are 
from the vertices and are equal to the angle deficits, the 
genus of the iso-density contour at 5c can be obtained 
by the surnmation of the angle deficit over all vertices 
(jGott et a]||l"98l iWeinberg et al.lll987D . In this way, the 
genus of the iso-density contours is measured indepen- 
dent of the two-point statistics such as the power spec- 
trum or correlation function, though mathematically we 
know that they should be related by Eq. ^ in the Gaus- 
sian case. 

As the structures grow, the fluctuations would be- 
come non-gaussian gradually, and the genus is no longer 
solely determined from the nonlinear power spectrum 
of the structure. However, the amplitude of the genus 
curve A does not change too much in the weakly non- 
linear regime, because the genus is a topological indica- 
tor which is independent of simple growth of clustering 
without merger. A general correction to Q{y) for the 
non-Gaussian field have been derived up to the secon d 
order perturbation (jMatsubara II1994I: lMatsubarall2003D 
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where iJ„(i^) = {-l)"'e''^^^{d/di')'^e-"^/'^ is the n-th or- 
der Hermit polynomials, S, T and U are the third order 
moments of 5 and its gradient VS. At the median density 
threshold = 0), these corrections are zero, so the am- 
plitude of the genus curve does not change. In practice, 
we calculate A from a set of points measured between 
V = —1 and +1. Since H2n+i{v) are all odd functions 
in A would not be affected by non-Gaussianity up to 
the second o rder. At smaller scales, further investigation 
(|Park II2011D shows that the directly measured genus am- 
plitude is better conserved than the shape of the power 
spectrum (Eq. [2]). The insensitivity of the genus to the 
nonlinear structure formation is not an artifact of the 
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smoothing procedure, but reflects the more fundamental 
nature of topological pr operty of the r andom field. 

As demonstrated by iPark fc Kim I ()2010[ ). when the 
correct redshift-distance relationship r(z) is adopted, 
one would be using the same comoving volume V{z) = 
{D\/ H){z) and smoothing scale Rg at different epochs. 
Here Da{z) is the angular diameter distance and H{z) 
is the Hubble expansion rate. Therefore the data at dif- 
ferent redshifts actually enclose statistically almost the 
same amount of structures, and the amplitude of genus 
curve A would be almost the same at different redshifts. 
This can also be seen from Eq. ([2]), which shows that A 
actually measure the slope of the power spectrum around 
the smoothing scale Rg . Since in the linear regime of GR 
only the growth rate of the structure, not the shape of it, 
evolves, A would be conserved. However, if an incorrect 
r{z) due to an incorrect cosmology is adopted, both V{z) 
and Rg would be mis-estimated. Since the topology of 
the structure is not scale-free, the genus enclosed in a 
wrongly sized volume and smoothed with a wrong scale 
would lead to deviation from the actual one, in this case 



AY{z,Rg^^ 



)R%^Y — ■^x{Rg,x)Rg^x^ 



(6) 



where Y is the adopted cosmology parameters while X is 
the true cosmology. The smoothing scale Rg for different 
cosmologies are related to each other by 



{Rg,x/Rg,Y)^ 



{D\/H)x/{D\/H)y, 



(7) 



Utilizing this effect, topology of the large scale structure 
can serve as a standard ruler in cosmology. 

Since the genus would be affected by modifying the 
density perturbation through Eq. ([2]), we now consider 
the evolution of density perturbations in modified grav- 
ity. A wide range of modified gravity theories which sat- 
isfy the basic requirement of being a metric theory where 
energy-momentum is covariantly conserved, can be stud- 
ied in the s o called parametrized post-Friedmann (PPF) 
framework ()Hu fc: Sawicki|[2"007fl . For these theories, on 
superhorizon scales structure evolution must be compati- 
ble with background evolution, on intermediate scales the 
theory behaves as a scalar-tensor theory with a modified 
Poisson equation, while on small scales, to pass strin- 
gent local tests, the additional scalar degree of freedom 
must be suppressed. The evolution of linear perturba- 
tions in theories which satisfy these conditions can be 
characterized by a few parameters. With the Newtonian 
gauge temporal and spatial curvature scalar perturbation 
^ and 3>, we introduce 
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In the absence of anisotropic stress, $ = —5*, and 5 = 
for the GR case. However, in modified gravity g may not 
be zero. From causality considerations, on superhorizon 
scales the evolution of metric perturbation must be de- 
termined entirely by the expansion rate and g to the first 
order of kn, with 



' H' H" ' 



limfc^^o(^$" 

= 0{k%) = i/cfcffKn. (9) 
where ' = d/dlna, Vm ~ kn is the velocity perturbation. 



kn — k/Ha. However, to the next order, we can then 
introduce to characterize the curvature perturbation 
on superhorizon scales. The subscript C = <I> — Vm/kn, 
which is the curvature perturbation in matter comoving 
gauge, indicates the fact that the LHS of Eq. ^ equals 
to C in the matter comoving gauge. On subhorizon scales, 
the lensing potential satisfies the modified Poisson equa- 
tion. 



4ttG 
1 + /g 



a^/9„A™(fc), 



(10) 



where Am is the fractional density perturbation, fc 
parametrizes modification to Newton constant. Equiva- 
lently, we can introduce a new scalar field if which char- 
acterise the extra source felt by nonrelativistic particles. 



kH{k) = 
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aV™A„(fc) + -fcV(fc), (11) 



where (p{k) satisfies the equation 
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Thus in the subhorizon regime, the scale-dependent de- 
viation from the GR is described by <f{k). In sum- 
mary, in the PPF parameterization, the modified gravity 
model is characterized by g{a^k), fQ{a), f(^{a) and the 
super horizon-subhorizon transition scale cr- The evolu- 
tion of these functions depends on the specific model. 
The fluctuation power spectrum for such theories can be 
calculated by using the public Boltzmann code camb with 
the P PF module Q (Fang ct.al.2008a.b. : .Lewis fc Bridii 
I2OOI ). 

3. MODELS 
As a first example we consider the f{TZ) theory, 

-rU + fin) 
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where TZ is the Ricci scalar and Cm is the matter La- 
grangian density. For a given expansion history H{a), 
e.g. the effective dark energy equation of state We — —1, 
in order to explain the late-time acceleration, the form of 
f{TZ) can be determined from a second order differential 
equation 

6 3 

where p is the total energy density, f-ji and fun a-re the 
first and second derivatives of f{TZ) with respect to TZ. 
A given expansion history permits a family of f{TZ) func- 
tions, the additional degree of freedom is usually char- 
acterized by Bq, the present day value of the function 
B{a), which is the square of the Compton scale given by 
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B{) is thus a model parameter of the fiJZ) theory, which is 
to be constrained by structure growth. Given H{a) and 

^ |http://camb.mfo/ppf/| 



B{a), the metric ratio at superhorizon scale gsH can be 
obtained by solving the differential equation 
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and utilizing the relation ' j / = ( -<& - i?<i>')/(l - B). 
As shown by IHu fc Sawicki' (2007), the scale dependent 
g{a, k) can be well fitted by the interpolation function 



g{a,k) 



gsH + ggsicgkH)""" 
1 + {cgknY^ 



(17) 
= 2, with 

other PPF parameters = —l/3g, fa = fn and cr = 1. 



where ggs = -1/3, Cg = 0.71B^/^ and 
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Figure 1. The genus per smoothing volume as a function of 
smoothing scale Rg at a = 1 for /(T?.) theories with different Bo. 

Given PPF parameters of f{TZ) theory above, espe- 
cially g{a,k), the matter power spectrum at relevant 
scales can be obtained with the help of Eq. (ITTl [T^ . 
In Fig. (HI, we plot the genus amplitude of the genus as 
a function of the smoothing scale Rg at a = 1 for several 
f{TZ) models characterised by difi'erent Bq values. For 
models with greater Bq value, the genus is larger. 

We may also see how the amplitude of the genus curve 
varies as a function of redshift. In Fig. ([2]), we plot 
the redshift evolution of A for a fixed smoothing scale 
Rg ~ 15Mpc//i. Compared with the horizontal line, 
which corresponds to GR in the bottom of the figure, 
the amplitude of the genus curve exhibits a strong time 
dependence for the f{Tl) theory with a total variation of 
about 10% for Bq > 0.01. 

From these figures it is obvious that for larger Bq, the 
deviation from GR is greater, i.e. the variation of the 
genus is stronger in the f{TZ) theory. By inspecting the 
perturbative variables in detail, with other PPF parame- 
ters we find that this is mainly due to the monotonically 
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Figure 2. The redshift evolution of genus amplitude for /(7^) 
theories v^ith different Bq. 

increasing deviation of $_ from the GR prediction with 
respect to the wavenumber k. The extra scalar degree 
of freedom ip — —f-ji enhances the gravitational force 
at smaller scales, and therefore increase the slope of the 
clustering. 

Next we consider the self-accelerating Dvali- 
Gabadadze-Porrati (DGP) braneworld model. In 
this model, our universe is a (3-1-1) dimensional brane 
embedded in an infinite Minkowski bulk, it can be 
described by the following action 



n 



(5) 
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n 
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167rG, 
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where 7?.(i), G'(i), i = 4, 5 is the Ricci scalar and New- 
ton constant in the b rane and the bulk respectively. At 
the background level, iDvali et al.l (|2000D showed that the 
consequent Friedmann equation 
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leads to a de-Sitter phase at late time characterized by 
the crossover scale Tc = G(5) 720(4). 

Above the horizon, the evolution of the perturbation 
can b e solved by iterative scaling method (| Sawicki et al.l 
l2007f ). and the metric ratio gsH is well fitted by function 



9SH{a) 



8Hrc - 1 
In the quasi-static regime. 



0.51 



Hrr - 1.08 



5Qs(a) 



1 H' 

1 - 2Hrc ( 1 + 



(20) 



(21) 



Thus the scale dependent g(a, k) can be described by the 
same interpolation equation Eq. (I17p with Cg = 0.4, Ug = 
3. Other PPF parameters are Jq = 0.4(7, Jg = 0, cp = 1- 
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From the calculation, we find that the major scale- 
dependent modification to the growth of the structure 
is in the superhorizon regime {k/aH <C 1). With a 
Gaussian smoothing scale of 15Mpc//i, only about 0.07% 
deviation of A arc found compared with GR, which 
can hardly be distinguished by the observation in the 
near future. In this case other techniques such as the 
ISW effect or the growth rate are i nore viable. How- 
ever it should also be not iced that ()Scoccimarrol 120091 : 
iChan fc Scoccimarrol l200"9l ) with the help of Vainshtein 
mechanism which bring the gravity back into GR at small 
scales, the nonlinear evolution can also induce significant 
scale-dependent deviations at relevant scales. We leave 
this topic to our further study. 

At last we consider more generic models in the context 
of PPF, this allows the method to be applied to other 
possible modified gravity models, and also the analysis 
would help us to see for whick kind of model the topo- 
logical method is more sensitive. We note that amongst 
the PPF variables, faia) and /({a) connect the metric 
fluctuations to the matter fields and are only functions 
of time. Rescaling the amplitude of the matter fluctu- 
ations without changing the shape (i.e. slope) of the 
power spectrum would not affect the genus, so here we 
do not need to consider them, as the topological mea- 
surements would be insensitive to these. Similarly, for 
the wavenumber range corresponding to smoothing scale 
of the order O{l0)Mpc/h relevant to large-scale genus 
topology study, the mode is in the sub-horizon regime, 
so the variable gsH which describes the modification to 
the superhorizon metric can also be excluded from our 
consideratiorL 

Following iHul (|2008[ ). we adopt the scale depen- 
dent g as an interpolation function between the super- 
horizon regime gsH and the quasi-static regime ggs (c.f. 
Eq. ([TTI) . Both Cg{a) and the power index Ug can affect 
the behavior of the transition between the two regimes. 

Eq. (|17p shows that for each Fourier mode the time- 
dependent factor Cgkn = Cgk/{aH) determines the mo- 
ment and scale at which the gravity deviate from the GR 
case. In this simple model, limfc^.o g{a, k) = 0, therefore 
(p would be modified when k ^ kt = {aH)/cg. 

First let us consider the case where Cg is a constant. 
The transition scale would first slowly decrease as kt ~ 
a~^/^ during the matter dominated era, then start to 
increase during the era of accelerated expansion. In order 
to get significant effect on scales of fc ~ C'(0.1)/i/Mpc or 
above, Cg should be at least larger than {aH)/k, which 
is about 0(0.01). 

As an example, in Fig. ([3|), we plot the genus per 
smoothing volume as a function of Rg for various ggs, 
while keeping gsn — 0, Cg{a) — 1, and Ug = 2. Com- 
pared with the GR case (i.e. ggs = 0), the difference 
is quite apparent. At any time, the metric deviation 
Eq. ([T7|) is a monotonic function of the wavenumber, 
this induce the same monotonicity in the deviation of 
from GR. Specifically, negative value of gqs and there- 
fore the positive extra force —V(p, will increase $_ as well 
as the slope of power spectrum at k > kt, and ultimately 
rise the value of the genus. On the other hand, positive 
9QS behaves in the opposite way. Nevertheless, owing to 
the slow evolution of kt, the variation of A between the 
present and z = 5 is less than 1% for \gQs\ < 0.5. 
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Figure 3. The genus per smoothing volume as a function of 
smoothing scale Rg at a = 1 for the phenomenological model of 
different ggg with Cg{a) = 1. 
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Figure 4. Deviation of (upper) and power spectrum P{k) 
(low er) compared wi th GR, assuming g(a, k) has the form of 
Eq. ^ and Eq. Cg{0) = 0.1. 



Next, we consider models in which Cg varies with time. 
Here the deviation scale kt of the theory can vary signif- 
icantly with time. As a toy model, let us assume that 
Cg (a) depends on the square of the scale factor 

Cgia)^CgiO)-a', (22) 

and parametrize the model by the value of Cg(0). In this 
case, kt will always decrease, and for a fixed wavenumber 
the absolute value of the deviation \g\ will always increase 
with time. Consequently, as illustrated in Fig. (|4]), $_(/c) 
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(upper panel) and the power spectrum (lower panel) will 
deviate from the GR case progressively. 
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Figure 5. The redshift evolution of gen us ampli tude for phe- 
nomenological models described by Eg. I ll 71 1 and Eq. ll22l l. The up- 
per panel shows models with different ggg while keeping Cg{0) = 
0.1. The lower panel shows models with different different Cg(0) 



while keeping ggg 
case in each panel. 



: 0.5. The dashed line corresponds to the GR 



In Fig. ([5]), we can see the evolution of the genus for 
different model parameters. In the upper panel, gqs is 
varied, taking values from -0.5 to -0.1 while Cg = 0.1 a^. 
As expected, for the cases with larger value of gqs, the 
genus amplitude take a faster raise. In the lower panel. 
We show the result of varying the value of Cg{0), while 
keeping ggs — 0.5 fixed. In this case, increase Cg(0) 
would move kt to larger scales, and shift the onset of 
the deviation to an earlier epoch, so the genus should 
deviate from GR further and earlier. However, when 
Cg(0) > C(l), the modification would move to superhori- 
zon regime, then it would become saturated with a nearly 
constant genus. 

4. PROSPECTS OF OBSERVATION 

The genus curve of iso-density contours could be mea- 
sured with future large scale structure surveys. In the 
following, we consider how such surveys could be used to 
constrain the modified gravity models with the topolog- 
ical measurements discussed above. 

The variance of the actual genus measurement is usu- 
ally estimated with the help of simulation. For the pur- 
pose of making forecast, an analytical calculation is more 
desirable. In principle, the uncertainty in the topologi- 
cal measurement is different from the uncertainty in the 
power spectrum of the large scale structure, this is espe- 
cially true for the general (non-Gaussian) case. However, 
as discussed previously, the Gaussian assumption should 
be a reasonable approximation for the purpose of making 
forecast, as far as the large scales are concerned. In the 
appendix, we estimate the minimal amount of variance 
cr^ (Eq. [28| by propagating from the uncertainty in the 



power spectrum tTp(fc). We also compared our analyt- 
ical es timation with the measurement from iGott et all 
(|2009f) . who utilized two volume-limited subsamples of 
luminous red galaxies (LRG) to measure the genus statis- 
tics: a dense shallow sample at 0.2 < z < 0.36 with 
smoothing length Rg — 21/i^^/Mpc, and a sparse deep 
sample at 0.2 < z < 0.44 with Rg = Sih^^/Mpc. Their 
result is consistent with Gaussian distribution with the 
amplitude A = 167 A ± 7.0 and A = 79.6 ± 6.0 respec- 
tively. By assuming a reasonable bias factor b ^ 2, our 
formula Ea. (|28l) gives aA/A ^ 4.5% for the shallow sam- 
ple and aA/A ~ 7.4% for the deep sample, which are 
very close to the measured values (4.1% and 7.5% re- 
spectively). This shows that Eg. ipSl) could indeed give 
good analytical estimate of a a ■ 

Throughout this section, we assume the Gaussian 
smoothing scale Rg = 15Mpc//i unless explicitly empha- 
sized otherwise. At high redshifts the non-linear scale 
decreases, one can use smaller Rg which makes the un- 
certainty of the genus much smaller. 

Taking O = In^ as the observable, the Fisher matrix 

is 



E 

redshift bins 



dpi 



do 



1 



CMB, ij, 



(23) 



where Pi include 51^/1^, r2;,/i^, /ig, n^, In^l^, r, We, as well 
as other modified gravity parameters, e.g. Bq for f{Tl) 
theory. The Planck prior is added effectively by the 
contribution i^cMB, which helps to break the degenera- 
cies between the various parameters. The summation is 
over different redshift bins with size of Az — 0.1. For 
cosmological parameters {flmh'^,ho,'We), both the per- 
turbations and background evolution contribute to the 
constraints. Therefore, derivatives in Eg. (1^5)) should be 
carefully calculated by taking into account the scaling 
relation 



^obs(^e) — Atruci^^^^ Rg) 

A(Z): 



(24) 
(25) 



Here DA,ref and iJref are the angular diameter distance 
and the Hubble expansion rate evaluated in the reference 
cosmology, which is used when reconstructing the posi- 
tion from redshift, and for simplicity we assumed it as 
the same as the fiducial cosmology. 

We consider the measurement of large scale struc- 
ture topology in extended redshift ranges with a num- 
ber of optical/near IR galaxy surveys at z < 2, 
such as the LAMOSl B (see als o [Wag gLet al.' (20091)), 
BOS^I, WFM OS(Bassett et al.l [2005I) . JDElvQ and 
BIGBOSS( Sch legel et al.l 120091) survev. as well as a few 

21cm intensity mapping experiments at < z < 5, 

e.g. the CRT dChang et al.l l2007t fSeo et al.l [2009h and 
the MWA0. The survey parameters we adopt are listed 
in Table ([1]). Needless to say, these parameters are only 
prelimnary estimates based on current planning, the pa- 
rameter for the actual projects are subject to change. 



^ http://www.lamost.org/ 

^ http:/ /cosmology. lbl.gov / BOSSA http://www.sdss3.org/ 
* http://www.jdem.gsfc.nasa.gov/ 
^ http://www.MWAtelescope.org/ 
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Table 1 

Parameters of future large scale structure surveys. 



Galaxies Survey 



redshift range 



sky area 



BOSS 
LAMOST 
WFMOS 

JDEM 
BIGBOSS 



0.2 < z < 0.7, 
0.0 < 2 < 0.7, 
0.5 < 2 < 1.3, 
0.5 < 2 < 2.0, 
0.2 < 2 < 2.0, 



10000 deg^ 
8000 deg2 
2000 deg2 
10000 deg2 
24000 deg2 



21cm Intensity Mapping 


redshift range 




other parameter 


integration time (hour) 


MWA 


3.5 < 2 < 5, 


Na 


= 500, As = Wtt'^ deg^ 


1000 


MWA5000 


3.5 < 2 < 5, 


Na 


= 5000, As = 167r2 deg^ 


4000 


CRT 


< 2 < 2.5, 


L 


= 100m, W = 15m X 7, 


10000 




10-» 10-» 0.0001 0.001 0.01 



Bo 

Figure 6. Likelihood distribution on Bq of /{TZ) theory for vari- 
ous optical/IR surveys. Two fiducial value (Sg = 0.01 and 0.0001) 
are shown in the figure. 



In Fig. we illustrate the likelihood distribution over 
Bq of the f{Tl) theory at fiducial values Bq = 10~^ and 
10^'' calculated for various galaxies surveys. As one of 
the projects under planning, the BIGBOSS can provide 
the most stringent constraints, the l-cr error on Bq is 
0.0039 at Bq = lO'^ and 5.38 x lO'^ at Bq = 10"''. Even 
for surveys which are on going or will start in the near 
future (e.g. BOSS and LAMOST), one could also gain 
considerably rich information about the redshift evolu- 
tion of the genus amplitude. The l-cr errors at the same 
fiducial values are 0.011 and 1.1 x lO"'' for LAMOST, 
0.010 and 1.04 x lO"'' for BOSS respectively. Notably, 
the constraining conclusion drawn from particular sur- 
vey depends on the true cosmology, not only quantita- 
tively but also qualitatively. For example, the WFMOS 
survey, which is deep but narrow and is more powerful 
than BOSS and LAMOST survey at Bq = 10~^ become 
insignificant at Bq = 10""*. This is because when Bq 
is small, the amplitude of genus at high redshift would 
quickly approach to the value of GR (c.f. Fig. [5]) and 
become nearly indistinguishable for a narrow survey. 

Unlike galaxies surveys, the 21cm intensity mapping 




10-» 10-' 0.0001 0.001 0.01 



^0 _ 

Figure 7. The same plot as Fig.(|6ll for 21cm surveys. 

experiments have relatively poor angular resolutions, for 
the CRT 

R,,s{z) = r{z)-^ (26) 

can reach to the order of 10 ^ 30Mpc//i at higher red- 
shift. Here, r{z) is the comoving angular diameter dis- 
tance at redshift z, A is the received wavelength and 
LcRT is the length of the cylinder of the GRT telescope. 
Therefore the smoothing scale Rg when measuring the 
genus should be larger than i?rcs- 

We plot in Fig. ^ the likelihood of Bq for the 21cm 
intensity mapping. For the CRT, two different choices of 
Rg are considered here. One is the constant Rg model 
by assuming Rg — 2 x max{i?rcs('Z)}, which equals to 
60Mpc//i in this case. The ilgure shows that the con- 
straint of this model is not very stringent for the f{Tl) 
model. Another choice is to smooth the data with vary- 
ing Rg at different redshifts. We divide the whole red- 
shift range (0 < z < 2.5) into 4 intervals, and assume 
the same factor of two relationship between the smooth- 
ing scale and the angular resolution within each redshift 
interval. In this case, the l-cr error equals to 0.0060 at 
Bq = 10-2 and 6.93 x 10"^ at Bq = 10"''. On the other 
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hand, the angular resolution of the MWA and MWA5000 
is sufficiently good so that we assume the same smooth- 
ing scale as the optical/IR case. However, for their high 
redshift coverage, the constraints are still not very strin- 
gent, especially at Bq = 10"'*. 

Numerous efforts have been made to constrain the 
f{TZ) model, with various observations such as the CMB 
anisotropics, supernovae, BAO distance, weak gravita- 
tional lensing, galaxies flow and clusters abundance. The 
currently strongest constraints combini ng all of the data 
gives Bq < 1.1 X 10"^ at 95% C.L. ljLombriser et al.l 
|201C1( ). The main constraining power comes from the 
low redshift cluster abundance data, while the integrated 
Sachs- Wolfe (ISW) effects from galaxies-CMB cross cor- 
relation also pro vides a moderate cons traint of Bq < 
0.42 at 95% C.L.l jLombriser et al.l 120101). On the other 
hand, the galaxies power spectrum data( Tcgmar k et al.l 
[20061) toge ther with CMBt.Spergel et . al 2006) and 
supcrnovae ()Astier et al. I 120051 ) , only places an upper 
bound on Bq of order unity (jSong et al.l l2007b[ ). Al- 
though theoretical calculation shows significantly en- 
hanced growth of large scale structure in the f{TZ) 
models(Song et al. 2007^, uncertainties in the nuisance 
parameters such as the galaxies bias b and nonliearity 
parameter Qni substantially weakened its constraining 
power. Assuming similar survey parameter as those in 
ISong et al.l(|2007bD . we find the l-cr uncertainty from 
genus measurement is around 0.2 at Bq = 10~^. This 
demonstrates how the topological measurement which is 
insensitive to nonlinearity and clustering bias could im- 
prove the large scale structure constraint on the f{TZ) 
model, though we should note that in real measurement 
the error might be larger than statistical estimates given 
here. 




■ .7 / / 

0.04 I " ' ' I ' ' ' ' I ' I ' ' I I ^ ' ' I I ' ' ' ' I ' ' ' 
-0.6 -0.55 -0.5 -0.45 -0.4 



Figure 8. l-cr constraint on ggs and Cg{0), assuming Cg(a) = 
Cs{0)-a2 



In Fig. ([U, we also plot the 2D constraint on ggs- 
Cg{0) for the phenomenological model considered in 



Eq. (I17l22p . Similar to the case of f{TZ) model, the BIG- 
BOSS provides a very stringent constraint, the l-cr errors 
are (Tg^^ = 0.036 and o-Cg{Q) — 0.014 respectively. For 
the CRT, CTggs = 0.042 and 0-^^(0) = 0.028 (assume the 
redshift- varying Rg). (Tqqs — 0.18 and <Jc (o) — 0.19 for 
the LAMOST, and (Jg^^ = 0.17 and 0-^(0) = 0.18 for the 
BOSS. 

5. CONCLUSION 

The topological indicators such as the genus of the iso- 
density contours provides an independent way for charac- 
terizing the large scale structure, complementary to the 
more often used two point statistics such as the corre- 
lation function and power spectrum. A significant ad- 
vantage of the the topological measurement is that it 
is less susceptible to the non-linear evolution and bias 
of the large scale structure, thus reducing the effects of 
these potential systematic errors, and may therefore be 
an more reliable way to extract the information on the 
large scale structure. 

In this paper, we studied the topology of the large scale 
structure as a measure of the scale-dependent expansion 
history of the universe in models of modified gravity the- 
ories. In the modified gravity theory models, the struc- 
ture growth can be scale-dependent even in the linear 
regime, the amplitude of the genus curve varies signifi- 
cantly with redshift, and hence it can be used as a new 
tool for distinguishing the models from the case of stan- 
dard gravity (i.e. GR). We illustrated this for the f{TV) 
theory, DGP braneworld model as well as a phenomeno- 
logical m odel parametrized wit h the PPF variables intro- 
duced bv iHu fc Sawickil (|2007| ). We find that the genus 
curves for these models are modified and evolve with red- 
shift due to the scale-dependent growth effect, hence the 
genus curve can be used as an observable to distinguish 
the modified gravity models from the general relativity 
theory. Finally, using the Fisher matrix formalism, we 
also made forecast on the sensitivity of this test with 
some current or future optical/IR and 21cm redshift sur- 
veys, showing that the method is a competitive way to 
test modified gravity models. 
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APPENDIX: STATISTICAL UNCERTAINTY OF 
THE GENUS MEASUREMENT 

Since the genus amplitude A can be expressed as a 
function of the power spectrum (Eq. [2]) , we may approx- 
imately estimate the statistical uncertainty a a by propa- 
gating the error from the uncertainty of power spectrum 
o-p(fc). 

Rewriting the observable 

O = \n{A) = (3/2)[ln(7W4) - ln(7W2)], 
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where is defined as 

Mn = f d^k k"-^Px{k)W{kRg) 



telescope is dominated by the sky: 



(27) 



Px{k) denotes the power spectrum under consideration, 
then the variance cr^ = Var(C') is simply 



(If 



Var(X2) Var(X4) Coy{M2,M4) 



Ml 



Ml 



M2M4 



Here Var(7W„) = Cov{Mrn,Mm), and Cov(A^,„, A^„) 
is propagated from aj, 

Cov{M,n,Mn) = {M„r,Mn) - {M,n){Mn) 

= f dk dfi k^"'+''-^'>W^{kRg) aljk). (29) 

''s J 

where is the cosine of the angle between k and the line 
of sight, Vs is the survey volume. We have assumed the 
covariance matrix of power spectrum is diagonal. 



Cov(P(fci), P(fc2)) = Soiki 



^^2(27r)2 2 



(30) 



We need to estimate the statistical uncertainty in 
power spectrum measurements for the different experi- 
ments. For galaxy surveys, the statistical uncertainty of 
Pg(fc) per Fourier mode includes both the cosmic variance 
and the shot noise due to the finite number of galaxies: 



c^Pg {k) 



P. 



(31) 



As is often done in such estimates (e.g. iSeo fc EisensteirTI 
([2003,)), we assume the density of the galaxy n{z) satisfies 
nP{k = 0.2) = 3. 

For the 21cm intensity mapping experiment, e.g. the 
cylindrical radio telescope (CRT) , in addition to the cos- 
mic variance and the shot noise due to the finite number 
of galaxies, there is also the noise due to the foreground 
and receiver (ISeo et all 120091) ■ 
with 

1 



P21cm{k)=pi PHl{k) + -] + 



fcs(57;ky + T,)A/ 



Vr 



with ps = fcsgTsigA/, g = 0.8 is the gain, tint is the 
integration time and Vr is the volume of a pixel, Tgig is 
the average brightness temperature which is estimated 
as 



(32) 



with a conservative assumption for the neu_tral hydrogen 
fraction XHiiz)flH,o — 0.00037. Tgky and Ta are average 
sky and antenna noise temperatures, which are assumed 
to equal lOK and 50K respectively. 

For the high redshift 21cm experiment such as the 
MWA, we assume that the system temperature of the 



sys 



250[(l + z)/7]^-''K 



(28) p(r) 



and the observation time is 

tk = {AeUr.t/X^Mk^) 

at kj_, where n{k±) is the number of baselines which 
observe the transverse component of the wavevector. 
This can be calculated from the array configuration, 
here we model the antennas distribution of MWA as 



' , with radius = 750m and a fiat core of ra- 
dius rc ~ 20m. For the hypothetical follow-up of MWA, 
denoted as MWA5000, we have assumed = 2km and 
rc = 80m. 

The results of our estimates are given in Sec. 4. 
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